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ABSTRACT 

Expressions  of  the  sway  and  yaw  damping  coefficients  are  obtained 
for  submerged  slender  bodies  of  revolution  moving  at  a  constant  forward 
speed  near  a  free  surface,  derived  from  the  Lagally  theorem  for  unsteady 
flow  as  extended  by  Cummins,  The  singularity  distributions  which  generate 
the  flow  consist  of  a  basic  doublet  distribution  identical  with  that  for 
a  deeply  submerged  body  and  a  higher  order  distribution  to  account  for 
the  effect  of  the  free  surface.  Numerical  values  were  obtained  for  a 
spheroid  and  the  results  were  compared  with  the  strip  theory  as  well  as 
the  case  of  zero  forward  speed. 


NOMENCLATURE 
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Kq  -  g/U2 

k:  -  u2/g 


1  -  L/2 

I  -  (M  ,  M  ,  M  ) 

x  y  z 

N 


N' 

z 


amplitude  of  radiated  waves 

semi-axes  of  elliptical  cylinder 

integrals  given  by  equations  (24)  and  (25) 

maximum  diameter  of  body 

energy  of  radiation 

mean  energy  of  radiation 

force  vector 

X  /? 

Froude  number  «  u/(gL)  ' 
submergence  of  body 
Green's  function 

potential  of  a  longitudinal  doublet,  see  equation  (6) 

potential  of  a  lateral  doublet,  see  equation  (10). 

gravitational  acceleration 

integrals  given  by  equations  (48)  and  (49) 

integral  given  by  equation  (A5) 

spherical  Bessel  function 


added  mass  coefficients 
length  of  body 

doublet  strength 
moment  vector 
yaw  damping  moment 
yaw  damping  coefficient 
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n,  n 
pl*  p2 

q  ■  (u,  v,  w) 
R(x) 
r 
r" 

S(x) 

S 

o 

t 

U 

v 

o 

W  -  ^(L/g)1/2 
x,  y,  2 
Y 
Y* 

P 
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P 
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inward  normal  of  surface 
see  equation  (9) 

induced  velocity  by  external  flow 

radius  of  body 

yaw  velocity 

position  vector 

cross-sectional  area 

maximum  cross-sectional  area 

time 

forward  speed 
sway  velocity 

cartesian  coordinates 
sway  damping  force 
sway  damping  coefficient 
see  equation  (M) 
see  equation  (9a) 
density  of  fluid 
velocity  potential 


u  frequency 

[  ]  terms  containing  cos  ut  in  [  ] 

superscript  indicates  the  order  of  small  quantity 


INTRODUCTION 


In  1958,  the  Davidson  Laboratory  of  the  Stevens  Institute  of 
Technology,  under  the  sponsorship  of  the  Bureau  of  Ships  Fundamental 
Hydromechanics  Research  Program,  initiated  a  project  to  study  the 
lateral  motions  of  surface  ships  with  three  degrees  of  freedom  in 
oblique  waves.  At  that  time,  the  three-dimensional  effect  and  the 
forward  speed  effect  on  the  damping  coefficients  had  been  detected 
both  theoretically  and  experimentally  by  Haskind  [l],  Hanaoka  [2], 
Havelock  [3,  4]  and  Golovato  [5,  6]  for  ships  in  heave  and  pitch 
motions.  Since  similar  investigation  on  the  lateral  damping  coeffi¬ 
cients  was  then  lacking,  the  present  work  was  undertaken  in  order 
that  proper  correction  of  damping  coefficients  can  be  made  in  the 
study  of  lateral  motions  of  surface  ships. 

The  project  advanced  very  satisfactorily.  By  the  early  summer 
of  1959,  all  expressions  of  the  sway  and  yaw  damping  coefficients  of 
submerged  spheroids  were  obtained,  being  ready  for  numerical  computation. 
Unfortunately,  difficulties  were  encountered  in  the  computer  program¬ 
ming  and  the  work  was  therefore  temporarily  suspended  as  both  authors 
departed  from  the  Davidson  Laboratory  in  the  fall  of  1959.  Resumption 
of  the  study  was  not  made  until  the  summer  of  1960  when  one  of  the 
authors  returned  to  the  Laboratory. 

Having  overcome  the  difficulties  which  arose  in  the  computer 
program,  the  authors  later  learned  of  J.  N.  Newman's  recent  work  on 
the  damping  coefficients  of  a  submerged  ellipsoid  [7].  Although 
Newman  has  treated  a  more  general  problem,  the  present  study,  being 
based  on  a  different  approach,  is  by  itself  interesting. 

In  the  evaluation  of  damping  coefficients  of  oscillating  ships. 


2  - 


two  methods  are  usually  employed,  being  based  on  the  consideration  of 
the  wave  energy  radiated  from  the  body,  or  on  the  direct  calculation 
of  the  force  and  moment  on  the  body* 

The  energy  method  was  developed  by  Havelock  [3,  8,  9]  for  two- 
dimensional  damping  coefficients  as  well  as  the  three-dimensional  case 
of  zero  forward  speed*  Generalization  of  the  method  to  the  case  of  a 
moving  ship  was  made  by  Newman  [10,  11]*  The  advantage  of  this  method 
is  that  the  damping  of  a  submerged  oscillating  body  can  be  evaluated 
from  the  energy  carried  away  by  waves  generated  by  the  singularity 
distribution  which  is  determined  by  the  motion  of  the  body  in  an  infi¬ 
nite  fluid*  The  correction  of  the  singularity  distribution  due  to  the 
effect  of  free  surface  can  be  neglected. 

Haskind  [1],  Havelock  [4]  and  Newman  [12]  evaluated  the  heave 
and  pitching  damping  coefficients  of  thin  ships  by  integrating  the  pressure 
over  the  body.  Kaplan  and  Hu  [13]  calculated  the  damping  force  and 
moment  of  submerged  spheroids  at  zero  forward  speed  by  utilizing  the 
Lagally  theorem  extended  by  Cummins  [14]  for  unsteady  flows.  Although 
it  is  seen  in  Reference  13  that  a  higher-order  solution  of  singularity 
distribution  is  necessary  for  the  evaluation  of  the  damping  from  the 
force  calculation,  the  final  expression  does  not  show  any  difference 
from  the  result  obtained  by  energy  method*  Rjrthermore,  in  the  study 
of  the  motion  of  the  body,  one  has  to  determine  not  only  the  damping 
force  and  moment  on  the  body  but  also  the  force  and  moment  due  to  the 
Inertia  effect*  Therefore,  it  seems  more  convenient  and  systematic 
to  evaluate  all  forces  and  moments  from  the  same  approach* 

In  the  present  report,  the  damping  coefficients  were  evaluated 
by  applying  the  extended  Lagally  theorem  [14]  which  enables  one  to 


determine  the  force  and  moment  on  the  body  from  the  singularity  distri¬ 
bution  instead  of  the  pressure  integration  over  the  surface.  This 
offers  the  advantage  of  reducing  surface  integrals  to  line  integrals 
in  the  present  study  of  slender  bodies. 

Numerical  calculation  of  the  damping  coefficients  of  submerged 
prolate  spheroid  were  carried  for  various  speeds  and  frequency  parameters. 
The  results  were  compared  with  the  case  of  zero  forward  speed  as  well  as 
that  obtained  from  the  strip  theory. 


SINGULARITY  DISTRIBUTIONS 


The  origin  of  the  cartesian  coordinate  system  moving  with  the 
body  is  taken  on  the  undisturbed  free  surface  with  the  x-y  plane 
horizontal  and  the  z-axis  vertically  upward.  The  slender  body  of 
revolution,  fully  submerged  at  a  depth  of  f  below  the  free  surface, 
is  moving  forward  along  the  x-direction  at  a  constant  speed  U  and 
oscillating  at  a  frequency  u  in  sway  or  yaw,  its  center  being  on 
the  negative  z-axis. 

The  velocity  potential  $  of  the  fluid  satisfies  the  Laplace 

equation 

A0  -  0  ,  (1) 

the  boundary  condition  on  the  free  surface  (z  “  0) 


and  the  boundary  condition  on  the  body  surface 


where  n  is  the  normal  of  the  surface. 

It  is  well  known  that  the  potential,  for  slender  bodies  of 
revolution  moving  in  an  infinite  fluid,  may  be  represented  by  doublet 
distributions  on  the  axis  of  the  body.  The  potential  due  to  the 
forward  motion  of  the  body  may  be  represented  by  a  doublet  distri¬ 
bution  in  the  x-direction 

Mx^  (x)  "  ~~  S(x)U 


f 


(4) 
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while  the  lateral  oscillation  of  the  body  may  be  represented  by  a 
periodic  doublet  distribution  (x)  cos  wt  in  the  ^-direction 

where 


1  +  1^ 

■  ■"4y  S(x)  vq  for  sway. 


(x)  -< 


(5) 


1  +  k» 
4tt 


S(x)  rx 


for  yaw. 


In  the  above  expressions,  vq  and  r  are  the  amplitudes  of  sway  and 
yaw  velocities,  respectively,  being  assumed  to  be  small  as  compared 
with  the  forward  speed)  S(x)  is  the  cross-sectional  area  of  the  body, 
be  ing  assumed  to  be  small  as  compared  with  the  square  of  the  body 
length,  the  superscript  therefore  indicates  the  order  of  the  small 
quantity)  the  subscript  of  doublet  distribution  represents  the  orien¬ 
tation  and  k^,  k ^  and  k’  are  the  added  mass  coefficients  for  the  cor¬ 
responding  modes  of  the  motion,  depending  on  the  slenderness  of  the  body. 

In  the  case  of  bodies  moving  under  a  free  surface,  waves  are 
generated  and  potential  must  also  satisfy  the  boundary  conditions  (2) 
on  the  free  surface.  Havelock  [15]  has  found  that  the  potential  of  a 
constant  doublet  in  the  x-direction  with  unit  strength  at  a  point 
(h,  o,  -f)  below  the  free  surface  is  given  by 


G„(x,y,zjh,o,-f )  -  (-iy - y)  (x  -  h) 


n 

7 


1  2 
* 


k  ek^z  ^sln  [k(x-h)  cos  9]  cos  (ky  sin  0)  sec  0dkd9 
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*/ 


K  (z-f )  sec  © 


2  3 

cos  [K  (x-h)  sec  9]  cos  (K  y  sin  9  sec  9)sec  9d9 
o  o 


(6) 


where  r^2  *  (x  -  h)2  +  y2  +  (z  +  f)2  , 


*2  -  (x  -  h)2  +  y2  ♦  (z  -  f)2  , 
Ko  =«  g/U2 


(7) 


and  the  asterisk  attached  to  the  upper  limit  of  integral  indicates  the 
principal  part. 

The  potential  of  a  periodic  doublet  with  an  amplitude  of  unit 
strength  in  the  y-direction  below  a  free  surface  may  be  found  from  the 
Green's  function  satisfying  the  boundary  condition  (2)  on  the  free 
surface  since  a  Green's  function  represents  the  potential  of  a  source 
and  the  potential  of  a  doublet  can  be  obtained  by  differentiating  that 
of  a  source  with  respect  to  the  position  coordinate. 

Green's  function  satisfying  the  boundary  condition  (2)  on  the 
free  surface  as  obtained  by  Peters  and  Stokers  [16]  is  expressed  as 
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where 
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sec2  0  [l  -  (1 


M  cos  0)1/2]  2 
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Arc  COS  -j-Stt 
4wU 


for  4<»>u/g  -  1 
for  4wU/g  >  1 


(9a) 


Equation  (8)  represents  the  potential  of  a  pulsating  source  moving 
at  a  constant  forward  speed  under  a  free  surface.  The  potential  of  a 
pulsating  doublet  in  the  y-direction  can  be  found  by  differentiating 
equation  (8)  with  respect  to  and  is  expressed  as 


G  (x,y,zjh,o,-f|t)  -  (-ij - ij)  y  cos  ut 

/  r  r 

1  2 
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/  P22  P2(z“f) 

-2K  /  - e  sin  [p_(x-h)  cos  0  -  ut]  sin  (p_y  sin©) 

Pi“P2  2  2 

^  •  sec2©  sin  0  d  0 


Tf  2  /  -  \ 

p.  P^z-f) 

+2K  /  - e  sin  [p^(x-h)  cos  0  -  wt]  sin  (p^y  sin  0) 


TT 

>/ 


Pi“P2 


TT 


•  sec  0  sin  0  d  0 


(10) 


the  first  term  containing  l/r^  in  both  equations  (6)  and  (10)  represents 
the  potential  of  a  doublet  in  an  infinite  fluid  while  the  remaining  terms 
represent  the  potential  due  to  the  Image  system  and  waves.  Because  of  the 
velocities  induced  on  the  body  by  these  remaining  terms,  the  boundary  con¬ 
dition  (3)  is  infringed,  and  additional  doublet  distributions  are  therefore 
needed. 

The  velocity  induced  on  the  body  axis  by  the  image  system  and  waves 
of  the  doublet  distribution  may  be  expressed  as 


where 


(u(2>.  v<2>,  w<2))  -  -  *(2) 

T 

J 

J*(2)  -/  Mx1)(h)  (Gx  ”  dh 

J  r. 


(11) 


(12) 


represents  the  potential  due  to  the  image  system  and  waves  of 
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being  the  half-length  of  the  body.  In  the  above  expression,  it  has  been 
assumed  that  the  submergence  of  the  body  is  sufficiently  large  so  that 
the  potential  due  to  the  image  system  and  waves  is  at  least  one  order 
higher  than  that  due  to  the  doublet  itself, 

(2) 

Equation  (6)  shows  that  Gx  is  an  even  function  of  y,  vv  ' 
therefore  vanishes  and  the  additional  doublet  distribution  on  the  body 
axis  to  correct  the  boundary  condition  distorted  by  may  be  ob¬ 

tained  as 


V(3)  -  o  ,  M^3)) 


[(1  +  1^)  2\  o  ,  (1  +  k2)  w^]. 


(13) 


(2) 

The  image  system  and  waves  of  the  doublet  distribution  in¬ 


duces  a  velocity 

(u<3>.  v<3).  „<3>)  -  -  V‘3> 


y»o,  z»  -f 


(14) 


on  the  body  axis,  where 
t 

.  f  (G - cos  ut)  dh 

w  '  ^  y 


(15) 


represents  the  potential  due  to  the  image  system  and  waves  of 
An  additional  doublet  distribution 


10 


on  the  body  axis  is  therefore  needed. 

Neglecting  the  higher  terms,  the  doublet  distribution  on  the  body 
axis  which  represents  the  body  to  the  fourth  order  may  be  written  as 


(1)  „  uO)  *  u(4>  m(2) 


M  -  (Mv  '  +  My'  +  M, 


cos  wt  + 


«<4>. 
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M 
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(17) 


The  corresponding  velocity  potential  may  be  written  as 

4  4 
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(18) 


where 


for  sway 
for  yaw 


(19) 


and  represents  the  potential  of  a  constant  doublet  in  the  z-direction 
with  unit  strength  at  a  point  (h,  o,  -f)  below  the  free  surface.  In  the 
present  study,  only  the  lateral  damping  force  and  moment  are  of  Interest, 
it  will  be  seen  later  that  the  doublet  distribution  in  the  z-direction 
makes  no  contribution  to  the  lateral  force  and  moment  and  the  explicit 
expression  of  Gz  is  therefore  unnecessary. 
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SWAY  DAMPING  FORCE 

The  Lagally  theoremf  which  allows  one  to  determine  the  force  and 
moment  on  a  body  without  carrying  out  the  pressure  integration  over  the 
body  surface,  has  been  extended  by  Cummins  [14]  to  unsteady  flows*  The 
extended  theorem  has  been  applied  successfully  to  the  evaluation  of 
stripwise  damping  coefficients  for  heave  and  pitch  of  a  spheroid  at  zero 
forward  speed  [13]  and  will  be  therefore  employed  in  the  present  study* 

The  extended  Lagally  theorem  [14]  for  the  force  may  be  expressed 
as 

7  -  -  4trp  ff  •  V  q-  4wp  ^  ,  (20) 

where  F  is  the  force  acting  on  the  doublet  M  and  q  is  the  velocity 
induced  by  the  external  flow  on  the  doublet.  In  evaluating  the  force  on 
a  body  generated  by  doublet  distributions,  the  mutual  interaction  only 
contributes  internal  forces  which  will  nullify  each  other,  therefore  the 
velocity  induced  by  other  doublets  in  the  body  should  be  omitted. 

Now,  the  velocity  induced  by  the  expernal  flow  may  be  expressed  as 

q  -  (-U,  v^,  o)  -  ?  -  V  +  .  .  .  .  (21) 

where  v^,  an(j|  tf(3)  are  given  by  equations  (19),  (12)  and  (15), 

respectively,  and  the  ellipsis  dots  ...  represent  higher  order  terms  of 

velocities  induced  by  the  image  system  and  waves  of  the  doublet  dlstri- 
-(4) 

bution  SP  '. 

In  the  evaluation  of  the  damping  force,  terms  containing  cos  ut 
are  of  interest.  Substituting  equations  (17)  and  (21)  into  equation  (20) 
and  keeping  terms  up  to  the  fourth  order,  one  obtains  the  stripwise  sway 
damping  force  as 


12  - 


§  -  4„p  «W  [i^i]  *  4np  CO.  „t  -  4irp  t^f-1 

oy  c  dy*  ( 


dM 


(4) 


(22) 


where  derivatives  of  potentials  are  evaluated  on  the  axis  of  body  and  the 

symbol  [  ]  indicates  that  only  the  terms  containing  cos  ut  in  [  ]  are 
c 

to  be  taken.  The  total  danping  force  can  be  obtained  by  integrating  the 
above  expression  over  the  length  of  body.  Assuming  that  the  body  is 
symmetric  about  its  midship  section,  i.e.,  the  cross-sectional  area  S(x) 
is  an  even  function  of  x  ,  one  finds  that  the  sway  damping  force  is 
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•  C^(x,h|  0,  2,  3)  dh  dx  . 
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m  _  ,  n  „ 

9  tan  9  e 
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secm  ©  sin11  ©  e  0  cos  [KQ(x-h)sec  ©]  d©. 
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YAW  DAMPING  MOMENT 


The  extended  Lagally  theorem  for  the  moment  is  expressed  as  [14] 

TT  -  -  4n  prx  (M  •  ?)  q"  +  4rrpqxM  +  J  p^rxn  dS  (26) 

where  NT  is  the  moment  on  the  doublet  M  about  the  origin,  r  the 
position  vector  of  the  location  of  the  doublet,  n  the  inward  normal  of 
the  surface  dS  and  the  integration  in  the  present  problem  is  to  be 
carried  out  over  the  peripheral  surface  of  the  strip  which  contains  the 
doublet. 

The  moment  contributed  by  the  first  two  terms  in  equation  (24) 
is  the  "Lagally"  moment  while  the  moment  contributed  by  the  surface 
integral  is  due  to  the  time  dependence  of  the  flow. 

Substituting  equations  (17 )  and  (21 )  into  the  first  two  terms  in 
equation  (26),  one  obtains  a  stripwise  yaw  damping  moment,  attributed 
to  the  Lagally  moment,  as 


dN 


dx 


~  -  4npxM^  +  4itpx  M^2 ^  cos  cot 

xdy  c  y  dy 

f<2) 


,2^(2) 


-4np  (U  +  ?C--)  M^2 ^  cos  <ot  -  4irpU  [M^4^] 
3x  y  y 


+4trp  M^  [rx  cos  cot 
x 


:W(3 )  /o\ 

+  — ]  +  4rrp  Mv  '  rx  cos  cot, 

3y  x 

c 


(27) 


where  derivatives  of  potentials  are  evaluated  on  the  axis  of  body. 

The  surface  integral  in  equation  (27)  may  be  evaluated  approxi¬ 
mately  for  the  present  problem.  By  using  the  method  developad  by 
Cummins  [17],  one  obtains  a  stripwise  yaw  damping  moment 
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dN 


2tf 


z2 


dz 


-  px  R(x)  J*  $  cos  0'  d0‘] 


(28) 


where  R(x)  is  the  radius  of  cross-section,  ^  the  potential  evaluated 
on  the  pheripheral  surface  of  the  strip  and  0’  the  angle  in  the  plane 
of  cross-section  measured  from  the  y-axis. 

In  examining  the  expression  of  potential  0  given  in  equation 
(18),  it  is  found  that 

J 

p  -  +  J  -Eg.  dh  (29) 

-J  rl 

is  the  only  part  in  potential  <ji  which  contributes  to  the  damping  moment. 
On  the  peripheral  surface  of  the  strip,  equation  (27)  may  be  written  as 


+  J 

-J 


R(h)  cos  9' 

[(x-h)2+R2(x)]3/^  ^  * 


(30) 


Since  the  integral  of  the  second  term  in  equation  (30)  peaks  sharply  in 
the  vicinity  of  the  point  x  *  h,  the  value  of  the  integral  is  mainly 
determined  in  the  range  of  the  peak.  Furthermore,  the  body  is  assumed 
to  be  slender,  i.e.,  R  <  <  I,  the  integral  in  equation  (30)  may  thus  be 
evaluated  approximately  in  the  same  manner  as  that  in  solving  the  integral 
equation  for  the  doublet  distribution  of  an  elongated  body  of  revolution 
[18]  and  equation  (30)  becomes 
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where  the  factor  2/ (l+k^ )  is  used  for  taking  into  account  of  the  effect 
of  finite  length  of  the  body. 

Substitution  of  equation  (16)  into  equation  (31)  yields 


<f>'  -  ^  -  R(x)  cos  0*  .  (32) 

The  potential  p^3^  varies  slowly  on  the  peripheral  surface  of 
the  strip  since  the  variations  of  y  and  z  are  small  in  comparison  to 
the  length  of  body.  One  can,  therefore,  expand  in  a  Taylor  series 

about  the  point  (x,  o,  -f)  as  Cummins  did  in  the  evaluation  of  the 
moment  on  a  slender  body  under  regular  waves  [17],  This  leads  to 


(x,  R  cos  S',  R  sin  0'  -f) 


S^3)(x,  o,  -f) 


+  R  cos  0* 


ajzK3) 

dy 


+R  sin  0' 


~dz  ' 


(33) 

where  higher-order  terms  have  been  neglected  and  the  derivatives  of  j^3^ 
are  evaluated  at  the  point  (x,  o,  -f). 

Combining  equations  (14),  (32)  and  (33),  one  finds  that 

0'  =  j^3^(x,  o,  -f)  -  2R  (x)  v^3^  cos  0'  -  R  (x)  w^3^  sin  0'. 


(34) 

It  follows  that  the  yaw  damping  moment  given  by  equation  (28)  may  be 
written  as 


2pxS(x)  [• 


dv 


(3) 


dt 


(35) 
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Summing  up  equations  (27)  and  (35)  and  integrating  the  result  over  the 
length  of  body,  one  obtains  the  total  yaw  damping  moment 
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where 


S^(x,hjm,n,q)  -  [  f  +(-)*  f  ]  sinm  9  tan11  9  e  * 

1-1  J  J 


Tf 

7 


-2p.f 


Tt 
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H  j 

- - sin  [p.(x-h)  cos  9]  d9, 

P1-P2  i 


(37) 


32(x,h#m,n)  - J 


TT 

7  -2K  f  sec"9 

secm  9  sin11  9  e  0  sin  [K  (x-h)  sec  9]  d9, 

o 


(38) 


and  and  C2  are  given  by  equations  (24)  and  (25). 


19 


DAMPING  COEFFICIENTS 

The  sway  and  yaw  damping  coefficients  Y*  and  N^  may  be  obtained 
from  the  sway  damping  force  and  the  yaw  damping  moment  by  the  relations 

Y  ■  Y*  v  cos  ut 

and 

N  «  N'  r  cos  ait. 
z  z 

In  the  case  of  a  prolate  spheroid  for  which 

S  (x)  -  So  [1  -  (j)2]  (40) 

where  Sq  is  the  maximum  cross-sectional  area  of  the  spheroid,  the 
integrations  with  respect  to  x  and  h  in  equations  (23)  and  (36)  can 
be  carried  out  explicitly.  Using  the  integral  expression  for  the  Bessel 
function 

(£)  n  *  7  1 

Jn  ♦  1  <z)  "  — - I7T  [  (1  '  t2)"  fllZt  dt*  (41) 

?  nl  v'  J 

-1 


>  (39) 


one  can  easily  show  that 
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and 
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16  j2(z) 


-1 
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where 


Jn  (z) 


(^)1/2  J.  -  1 


n  +  2 


(z) 


(45) 


is  the  spherical  Bessel  function. 

Substituting  equations  (23)  and  (36)  into  equation  (39)  and  em¬ 
ploying  the  relations  given  by  equations  (42)  through  (44),  one  obtains 
the  sway  and  yaw  damping  coefficients  for  a  prolate  spheroid 
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The  singular  point  cos  9  =*  g/4<»>U  where  p^  =  p2  in  the  expression 

of  ^(m.n.qjr.s)  can  be  removed  by  an  appropriate  transformation  of  the 

variable  given  in  Appendix  I.  Numerical  evaluation  of  1^  and  I2  was 

* 

then  carried  out  on  digital  computers  for  a  spheroid  with  a  length-diameter 
ratio  *  8  and  submerged  at  a  depth  of  two-tenths  of  the  length  below  the 
free  surface.  The  resulting  damping  coefficients  were  plotted  in  Figures  1 
and  2,  together  with  their  values  calculated  for  zero  forward  speed  as  well 
as  those  obtained  from  the  strip  theory  (Appendices  2  and  3). 


♦ 

I2  was  evaluated  on  IBM  704  by  the  Applied  Mathematical  Laboratory,  David 
Taylor  Model  Basin,  Washington,  D.  C. 


DISCUSSION  OF  THE  RESULTS  AND  CONCLUSIONS 


The  lateral  damping  coefficients  evaluated  from  equations  (46)  and 
(47)  for  a  spheroid  as  shown  in  Figures  1  and  2  do  not  Indicate  any  singular 
behavior  at  the  point  uU/g  =  l/4  which  is  contrary  to  the  case  of  longi¬ 
tudinal  damping  coefficients  [4,  7,  11].  This  basic  difference  between 
the  lateral  modes  of  oscillation  and  the  longitudinal  modes  of  oscillation 
can  also  be  seen  from  the  Green's  function  given  by  equation  (8)  which 
represents  the  potential  of  a  pulsating  source  moving  at  a  forward  speed 
under  a  free  surface.  The  wave  terms  in  equation  (8)  contain  a  denominator 
p^  -  which  equals  sec^  0  (l  -  cos  0)  at  the  point  wU/g  *  l/4, 

therefore  the  integrands  behave  like  l/O  near  the  point  9=0  and  the 

integrals  are  singular. 

The  potential  of  a  doublet  in  the  longitudinal  direction  (x-direction) 
or  vertical  direction  (z-direction) ,  obtained  by  differentiating  that  of  a 
source  given  by  equation  (8)  with  respect  to  £  or  £  ,  has  the  same 

singular  point.  However  the  potential  of  a  lateral  doublet  (in  the  y-direction), 

obtained  by  differentiating  equation  (8)  with  respect  to  N  ,  brings  out  an 
additional  term  sin  9  to  the  integral  of  the  wave  terms,  which  then  does  not 
have  any  singular  behavior.  Since  the  flow  due  to  longitudinal  modes  of 
oscillation  can  be  represented  by  distributions  of  longitudinal  and  vertical 
doublets  while  lateral  modes  of  oscillation  can  be  represented  by  distri¬ 
butions  of  lateral  doublets,  the  former  thus  has  a  singular  point  at 
uU/g  =  l/4  and  the  latter  does  not. 

Physically,  longitudinal  modes  of  oscillation  generate  waves  which 
are  symmetric  about  the  median  plane  of  ship  and  the  amplitude  of  waves 
becomes  infinite  at  the  critcal  point  uU/g  =  l/4.  Havelock  [4]  inferred 
that  this  behavior  probably  arises  from  the  phenomenon  of  resonance.  In 
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the  case  of  lateral  modes  of  oscillation,  however,  waves  are  anti-symmetric 
about  the  median  plane  of  ship,  i.e.,  waves  generated  by  one  side  of  the 
ship  are  180°  out  of  phase  with  waves  generated  by  the  other  side  of  the 
ship.  It  is  therefore  believed  that  this  anti-symmetry  of  waves  smooths 
out  the  singularity  at  the  critical  point. 

Figures  1  and  2  also  show  that  the  three-dimensional  effect  and  the 
forward  speed  effect  are  rather  important  on  the  damping  coefficients.  In 
order  to  present  the  results  in  a  more  practical  form,  the  authors  have 
tried  to  obtain  three-dimensional  correction  factors  for  damping  coefficients 
at  various  forward  speed,  but  have  found  that  these  correction  factors,  un¬ 
like  the  case  of  zero  forward  speed  treated  by  Havelock  [3],  Vossers  [19], 
and  Kaplan  and  Hu  [ 13] ,  are  highly  irregular  with  regard  to  the  variation 
of  frequency  parameters  and  it  is,  therefore,  impossible  to  draw  any  useful 
conclusions . 

The  behavior  of  the  damping  coefficients  in  the  present  study  agrees 
very  well  with  that  in  Newman's  work  on  the  ellipsoid  [7]  except  for  low 
frequencies  at  high  Froude  numbers.  The  discrepancies  may  have  resulted 
from  the  fact  that  the  boundary  condition  (3)  on  the  body  surface  is  satis¬ 
fied  on  the  mean  position  of  the  surface  in  the  present  study,  but  it  is 
satisfied  on  the  instantaneous  position  of  the  surface  in  Newman's  work. 

As  Newman  mentioned  in  his  oaper,  to  satisfy  the  boundary  condition  on  the 
instantaneous  position  of  the  surface  is  equivalent  to  including  the  "effect 
of  angle-of-attack"  due  to  forward  motion  and  the  effect  due  to  singularity 
distributions  located  on  an  oscillating  surface.  The  deviation  between 
Newman's  work  and  present  result,  therefore,  shows  that  these  effects 
should  be  taken  into  account  for  cases  of  low  frequencies  at  high  Froude 
numbers.  However,  the  question  arises  whether  the  other  "effect  of  angle- 
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of-attack"  due  to  forward  motion  that  generates  the  "lift"  on  the  body  should 
also  be  included  in  the  study  of  damping  coefficients.  It  is  well  known 
that  a  body  placed  at  an  angle  of  attack  in  a  stream  exerts  a  lift  force 
which,  in  the  idealized  theory  of  a  perfect  fluid,  is  generated  by  vortices 
shedding  away  from  the  body.  If  the  body  is  also  oscillating,  the  unsteady 
lift  contains  a  damping  component  resulting  from  the  energy  shedding  away  in 
the  oscillatory  vortex  wake.  Since  this  oscillatory  wake  also  generates  waves 
if  the  body  is  moving  under  a  free  surface,  it  is  therefore  important  to  find 
out  whether  this  has  any  appreciable  effect  on  the  damping  coefficients.  One 
of  the  authors  is  now  studying  the  case  of  a  thin  ship  which,  fortunately,  can 
be  treated  by  the  existing  lifting  surface  theory.  The  result  will  be  pre¬ 
sented  in  a  later  report. 
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APPENDIX  I 


A-l 


Evaluation  of  the  Integral  I^(m,n,q,r,s)  in  the  Text 

Since  the  integrand  of  I1(m,n,q,r,s)  given  by  Equation  (48)  has 
a  singularity  at  the  point  cos  9  =  g/4<oU  where  -  P2  vanishes,  an 
appropriate  transformation  of  the  integral  is,  therefore,  necessary  in 
order  that  a  mechanical  quadrature  can  be  performed. 


Defining  a  Froude  number 

F  =  u/CgL)1/2 

(Al) 

and  a 

frequency  parameter 

W  =  w(-)1//2, 
g 

(A2) 

where 

L  =  21  is  the  length  of  body. 

one  may  then  rewrite  equation  (9) 

as 

pik 

/■  “  (l  -  p)  sec  8  , 

p2L  4 

(A3) 

where 

p  -  (1  -  4Ftf  cos  9)1/2. 

(A4) 

Since 

the  difference  between  p^  and 

p^  occurs  only  in  the  sign  attached 

to  p,  it  is  sufficient  to  deal  with  the  integral 


it  q 

?  P1 

I' (m,2,q,r,s)  *  f  sec™  9  tan  9  - 

1  J  P l“P2 

TT 

2 

—2p  f 

•e  j*  (p^l  cos  9)  j®  (p^  cos  9)  d0  ,  (A5) 

wnich  may  be  written  as 


A-2 


I*(m,2,q,r,s)  •  L1_q J 


sec  q  tan2  e  Si* til  [-m£  ^ 
4 P  4F2  cos2  9 


m 


7 


„  U+P)  .  f  sec2@ 

•  e  2F2  T  C(-^)2  1 


r) 


8  cos  9J 

•  %  1  <«> 


The  singular  point  at  p  =  0,  i.e.,  at  cos  9  -  g/4uU  ,  can  be  removed 
if  one  chooses  p  as  the  variable  of  intearation.  It  follows  that 

(l+AFU)1^ 

Ij(m,2,q,r,s)  -  22m  +  2q  “  1  -L  1'q.Fn,+1.wm+2q“1  /*  ( 1- p )“2q ( !- p2 )“m 


1  /4W  v2  f 
VTTBV  *  T 


ri  ,1-p  ^2-,l/2  "  7  *  L  .r  r  W(l+p)n  ,s  r  W(1+P)t  JO 

*  U  -  (4F§~)  ]  •  -Jj.  [2?CTy]  ^2  dP 


(A7) 

where  the  integrand  has  no  singularity  in  the  range  of  integration. 

Integral  containing  pq  in  equation  (48)  can  be  treated  analogously. 
Equation  (48)  can  finally  be  written  as 


I1(m.2,q,r,s)  -  22ni+2q~1  .  L1_q  .  Fm+1  .  W1"*2^1  [  J 


-n 


-I 


(1+4FW)1/2 


-(1+4FW)1/2 


/ 


]  o-p)"2q  o-p2)""  n-(ii!if)]l/2 


1  (Sr)24 


"  7  L  .r  r  W(l+P)n  ,s  r  W(l+P)n  Jn 

•  •  •  h  trFtr=fr3  h  liFn^r]  dp  * 


(AS ) 


A- 3 


where 


for  4FW  >  1 

(A9) 

for  4FW  -  1 


Using  a  ten-point  Legendre-Gauss  quadrature  [19],  the  above 
integral  was  programmed  for  the  digital  computer  IBM  650  with  f/L  •  0.2. 
The  results  were  then  used  for  the  evaluation  of  the  damping  coefficients 
given  by  equations  (46)  and  (47), 


APPENDIX  2 
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The  Case  of  Zero  Forward  Speed 


If  the  forward  speed  of  the  oscillating  spheroid  is  zero,  the 
velocity  potential  for  the  motions  of  the  fluid  satisfies  the  boundary 
condition 


(AlO) 


on  the  free  surface  z  *  0.  The  potential  of  a  unit  doublet  satisfying 
the  above  boundary  conditions  may  be  obtained  by  differentiating  the 
Green's  function  [ 16]  with  respect  to  the  position  coordinate.  It 
follows  that  the  potential  pf1  may  be  written  as 


J  (ka«)  K.(z-f)  ,_v 

-  e  1  k  dk]  M'2' (h)dh 

k  -  K  y 


+  i  e 


iut 


2rr 


“J 


y  e 


K:(z-f) 


I 


V^a*) 


M^2)(h)dh 


(All) 


where 


o'2  -  (x-h)2  +  y2 .  (A12) 

2 

®  /g.  (A13) 

is  the  Bessel  function  of  the  first  order  and  is  the  doublet 

strength  given  by  equation  (5). 

At  large  distances  from  the  body,  terms  in  the  brackets  of  equation 
(All)  diminish  more  rapidly  than  the  last  terms.  The  real  part  of  the 
asymtotic  velocity  potential  may  be  found  as 


A-5 


.  i 

-  -2vKl2  e  1  J  sin  (wt  "  JCj_  a1  +  7  *)  My^OO  dh 

-t 

(A14) 

as  a’  — •>  00  • 

The  energy  radiated  from  the  body  may  be  determined  by  integrating 
the  product  of  the  pressure  and  the  radial  velocity  over  the  surface  of  an 
infinitely  long  vertical  cylinder  of  large  radius,  i.e.,  the  energy  of 
radiation  may  be  expressed  as 


(A15) 


where  o»  0  and  z  are  the  cylindrical  coordinates  with  a  the  radial 
axis,  9  the  angle  and  z  the  vertical  axis.  Since 


x  *»  a  cos  0  and 

one  thus  has 

a1  a  [(x-h)2  +  y2]1//2  =  a  -  h  cos  9 

and 

1  ~  1 

r’”  a 

as  a  -*■  CD  •  It  follows  that 


y  =  a  sin  9  , 

(Aid) 


(A17) 


(A18) 
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asymptotically,  and 
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E  «  lim  4wpK^  03  e  I  sin  0  [  I  cos  (o3t-K^a**K^h  cos  0  +  —  rr) 

a-*  as  ^  • 


2ir 


~l  .M^2\h)dh]2d0  . 

y  (A20) 


The  mean  energy  radiation  is  thus  found  to  be 

2tt 


_  3 

E  ■  2tt  pK^ '  os  e 


■2K,  f 


1  J*  (P2  +  Q2)  sin2  0  d  0 


(A21) 


where 


P  +  iQ  =■  J  M^(h)  e  1 


iK, h  cos  0 


dh. 


(A22) 


The  relation  between  the  damping  coefficient  and  the  mean  energy  radiation 
is  given  by 


E  -■ 


1  V.  2 

j  Y'  v 


1  2 

2*1  r 


for  sway 


for  yaw. 


(A23) 


one  finds  that  the  damping  coefficient  at  zero  forward  speed 


-2K.  f 


2tt 


'  -  ^  PK:3  03  (l+k2)2  e  1  J  (P:2  +  Qx2)  Sin2  0  d  0 


(A24) 


and 


-2K,  f 


2tt 


^z  =  pKl3  “  (1+k')2  e  1  J  (P22  +  Q22)  sin2  0  d  0 


(A25) 


where 
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-J 
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iK.h  cos  9 

S(h)  e  1  d  9 


*  «2  "/ 


iK.h  cos  9 

hS(h)  e  d  9. 


-J 


(A26) 


(A27) 


For  a  prolate  spheroid,  equations  (40)  and  (41)  may  be  employed  and  the 
above  expressions  become 


P1  +  iQl  “ 


4  (J) 


1/2 


S  J 

o 


(K^l  cos  9) 


§7?  J5/2  V  “S 


(*28) 


and 

4  (?)  l^2  */ 

P2  *  1Q2  '  1  (V  ,'c'  V*  (V  COa  S)-  (A29) 

The  damping  coefficients  given  by  equations  (A24)  and  (A25)  were 
evaluated  numerically  for  the  same  spheroid  treated  in  the  text,  and  the 
results  were  plotted  in  Figures  1  and  2, 
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Strip  Theory 

Havelock  [8]  has  found  that  the  wave  amplitude  at  infinity,  generated 
by  an  infinite  elliptic  cylinder  oscillating  horizontally  at  a  depth  of 
below  the  free  surface,  is  given  by 


A  ■  2"Kl»l  IT  (^)1/2  e_Klf  •  [K1  <*12  -  *22>  1/21- 


(A30) 


where  a^  and  a^  (a^  >  a^)  are  the  semi-axes  of  the  cylinder. 

For  a  circular  cylinder,  the  wave  amplitude  may  be  obtained  by  ex¬ 
panding  the  modified  Bessel  function  I  in  equation  (A30)  for  small 
2  2 

~  a2  and  setting  a^  =  a^  “  R  after  the  cancellation  of  a^  -  a^ 
in  the  denominator.  This  yields 


0  ,  v  -K.f 
A  *  2irK.  V  -1  e  1  . 

1  to 


(A31) 


The  mean  energy  radiation  per  unit  length  of  the  circular  cylinder 
2 

is  p  g  (0  A  /2K^  .  Assuming  that  this  holds  for  each  strip  of  the  spheroid, 
then  the  mean  energy  radiation  per  unit  length  of  a  strip  dh  is 

_2K  f 

dT  =*  2  p  u  vq2  e  1  S2(h)  dh  ,  (A32) 

for  sway  oscillation  and 

-2K  f 

dE"  **  2  p  K^2  to  r2  e  ^  h2  S2(h)  dh  (A33) 

for  yaw  oscillation.  Integrating  equations  (A32)  and  (A33)  over  the  length 
of  the  body  and  substituting  equation  (40),  one  obtains  the  damping  co- 
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(FIGURE  I.  SWAY  DAMPING  COEFFICIENT  FOR  THE  SPHEROID 

L/D  *  8,  L/f  *  5 
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FIGURE  2.  YAW  DAMPING  COEFFICIENT  FOR  THE  SPHEROID 
L/D  *  8, L/f  *  5 
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